





Постановка задачи. В работе | | изучена возможность использования фильтрации с помощью 
аналитического продолжения в верхнее полупространство в качестве способа разделения воз- 
мущающих источников на различньх глубинах. В зтом исследований сделан ряд допущений 
об упрощенном рассмотрениий геологической средь, но главное, что следует держать в виду, 
что грави- и магнитнье данньге не содержат явно информации о глубине источников аномалий. 
Вследствие решения соответствующей прямой задачи всегда можно создать любую модель на- 
блюденного поля по распределению значений намагничения (или плотности) на поверхности 
измерений. Но такое формальное математическое решения лишь в исключительньх случаях 
будет иметь геологически значимое содержание. 


Исходя из зтого в работе |1| сделано предположение о том, что наблюденное магнитное поле 
представлено в виде аддитивной модели: сумма полей регионального фона, остаточньтх анома- 
лий и агрегированного шума. При зтом, в качестве модели средьг использовано семейство взай- 
мно некоррелируемьгх распределений тонкослоистьх источников с целью решения прямой за- 
дачи - вьтчисления модели поля (получить зквивалентньй набор шаблонов полезного сигнала). 


Такой прим означаєт, что поля, вьічисленньге для разньїх значений глубин источников, имеют 
нулевую корреляцию. Ото условие отнюдь не вьшлядит таким ограничивающим, как может 
показаться на первьій взгляд, поскольку множество распределений неоднородности во вмеща- 
ющей среде, вообще говоря, могут иметь как положительную, так и отрицательную корреля- 
цию. Вследствие зтого допущение о нулевой корреляции рассматриваєтся всего лишь как част- 
ньшй случай нейтрального распределения. В физическом смьгсле, подобнье распределения пре- 
дставляют собой слои с зквивалентньїюми источниками, и к их обработке можно привлечь соот- 
ветствующие методьт фильтрации по зквивалентньєм источникам. 


Метод исследований. Обобщая вьшшеупомянутую работу ||, можно предложить следующий 
алгоритм восстановления потенциального поля по зквивалентному распределению источни- 
КОв: 

- | Для извлечения (отделения) поля регионального фона от источников, находящихся на глу- 
бине йо, нужно аналитически продолжить наблюденное поле на вьгсоту 2/0 по известньм со- 
отношениям. В волновой области оптимальньй фильтр для такой операции следующий: 
Е() з е 7"), Можно затем извлечь зто региональное поле из полного магнитного поля, 
чтобь получить разностное поле остаточньїх аномалий. В настоящее время имеєется доста- 
точно программньх средств, чтобьт осуществить подобную операцию с надлежаще подгото- 
вленньтми данньхми. 

- Для вьтчисления распределения намагничения в заданной области от слоя зквивалентньтх 
источников, расположенного между глубинами й; и й», нужно вьгтчислить разницу операто- 
ров аналитического продолжения вверх на вьгтсотах 2Л» и 2Лі, т.е., найти разность от пере- 
счета вверх сдвоенньх значений соответствующих глубин. В волновой области оптималь- 
ньй фильтр для такой операции: Б,(0) се )- є, 


Аномалий, имеющиєе источник внутри слоя зквивалентньїх источников, обьтічно вьіделяются 
более рельефно после фильтрации. Впрочем, в слабоградиентном мозайчном поле, осложнан- 
ном резкими структурньюми несогласиями, такая определінность теряется. Мьг рекомендуєм 
использовать альтернативньве методьт фильтрации (например, гармонический анализ |2|, либо 
фильтрьт Винера (31). 





В сообщении рассмотрен изложенньй вьше способ качественной фильтрации магнитньх ано- 
малий в приложений к обработке набора данньїх полного магнитного поля, полученньтх при 
детальном исследований археологического обьеекта в Ріпеуіеуу парке (исходнь6е даннь6е полу- 
ченьгт из репозитория у/ууму пес. поаа.соу). Здесь реализован позтапньгшй пересчет аномалий пол- 
ного магнитного поля вверх на вьгсоту 0.5, 1.5, и 4 м, соответственно (рис. 1) на равнинном 
участке территории парка 30х50 м (горизонтальньй масштаб дан в м, магнитная индукция - в 
нТ, контурьгт изолиний проведень через 25 нТ). Пересчет осуществляется с целью подавить ко- 
ротковолновой шум и вьщелить локальньг6"е аномалий. 
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Рис. І. Исходнье даннье (аномалия полного магнитного поля АТ в нтТ, верхний левьшй рисунок) 
и результать его продолжения вверх (слева направо), на 0.5 м; на І.5 м; на 2 м; на 4м. 


Чернье точки на рис. 1.1 обозначают расположение пунктов измерений с магнитометром. При 

зтом на рис. 2-5 нижний край участка обрезан из-за вьшоко-градиентного шума (близость к 

металлической ограде, исключено из анализа), а также, исключеньг (пропущенньєе точки) неко- 
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торьте наблюдения про профилях из-за низкого качества. Повторить наблюдения не представ- 
ляется возможньтм из-за изменения технологических условий (раскопь, обводнение и т.п.). Да- 
же в период измерений многие аномалиий из интересующего нас метрового диапазона, попали в 
проблемнь6е по условиям измерений зоньї. 


Магнитньєе аномалий, рассчитанньве от слоя зквивалентньх источников (рис. 2) в программе 
УїєпРгос 1.56 (С.В... Соорег, 2000) путем дифференцирования продленного на соответствую- 
щий уровень магнитного поля (рис. 1.2-5) между уровнями 0.5--1.5 м, 1-2 м, и 2-4 м, соответ- 
ственно. Зти слой заданньх источников априори содержат коррелируемье решения. Краткий 
качественньй анализ полученного поля позволяет предположить, что найбольший интерес с 
точки зрения планируемьх раскопок представляют аномалий, с центрами в координатах (по се- 
ти участка) х-8, у-20; х-з, у-29; х-30, у-22. Оти З аномалиий либо содержат заметное количест- 
во магнитного материала возле поверхности (2-4 м), либо вьжтзвань возмущающими телами, ра- 
сположенньми на значительной глубине (17-20 м под поверхностьк). 
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Рис. 2. Магнитнье аномалий, рассчитаннье от слоя зквивалентньх источников в слоях 
(0.5-1.5 м), (1-2 м), (2-4 м). 





Анализ результатов. Качественньй анализ исходного поля (рис. 1.1) и его трансформаций 
(рис. 1.2-5) подтвердил, что, в общем случае, продолжение вверх на величину, равную шагу по 
профилю (0.5 м), не даєт значительную "прибавку" в разрешений, но иногда значительно сни- 
жаєт вьісокочастотньтшй (коротковолновой) шум. 





Наложение сети точек наблюдений (рис. 1.1) на те же точки значений магнитного поля (рис. 
2.3), полеченньве при расчете прямого зффекта от слоя зквивалентньх источников на интерва- 
ле глубин 1-2 м (рис. 3.1) показало, что следует с крайней осторожностью трактовать 2 боль- 
шие аномалий на западном краю участка карть. Хотя большинство точек наблюдений и подт- 
верждают зти аномалий, их качественньй анализ показьгваєт большую вероятность глубинно- 
го происхождения аномальньх источников. К тому же, они сильно коррелируются с профиля- 
ми наблюдений, на которьїжх исключено много (1590) наблюдений с повьшшенньтм шумом. 


Последующее наложение значений магнитного поля, рассчитанньх на том же интервале глу- 
бин 1-2 м расположения слоя зквивалентньх источников, но с большей детальностью (рис. 
3.2, контур через 2 нТ) подтверждает предьщдущий вьгвод. Здесь точки указьтвают максималь- 
ньже значения горизонтального градиента поля. По их расположению можно качественно оце- 
нить положение границ возмущающих тел ниже поверхности наблюдений. В данном случає, 
особенности поля возле северного края картьг, менее вьтражень со сжатьтм контуром и подчер- 
кнутьт значительной величиной горизонтального градиента магнитного поля. Зто означаєт ме- 
ньшую вероятность содержательного истолкования (сомнительнь6е источники). 
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Рис. 3. І. Наложение сети наблюдений исходного магнитного поля на результат его трансфо- 
рмациий в поле зквивалентньх источников; 2. поле зквивалентньх источников с повьшенной 
детальностью. 


Альтернативньг-е исследования. В работе |4| исследована возможность восстановления моде- 
ли магнитньх свойств ниже поверхности наблюдений на основе решения обратной задачи маг- 
нитометриий с помощью адаптивного метода Кочнева (в программе АРСМ-3Д). Там успешно 
восстановлена послойная магнитная восприиймчивость для модели из 3 слобв. Оту методику 
можно рекомендовать для параллельного подечета и контроля вьічислений. 


Детальное исследование (5| обосновьгтваєт методику восстановления магнитньх свойств при- 
поверхностньх обьектов геологической средьт по данньм магнитометрии (модуль магнитной 
индукции В геомагнитного поля). В зтой работе исследованьт различньвев методьт разделения 
аномалий (усреднение; вьчисление вариаций; исключение регионального фона; аналитическоє 
продолжение вверх и вьтчисление остаточньх аномалий; фильтрация регионального фона). Ко- 
личественная характеристика зтих методов локализации магнитньх аномалий позволяєет уста- 
новить, что найболее предпочтительньтм является корреляционньєй метод вьшщеления локаль- 
ньїт аномалий, основанньй на тренд-анализе исходньгх данньжх. Он найменее подвержен пог- 
решностям, и рекомендован в |5| как лучшее средство разделения локальньгх и региональньтїх 
аномалий при обработке материалов магнитометрических сьбмок археологических обьеектов. 





Для оценки достоверности вьжявленньх аномалий использован приєм, описанньгй в |5|: форму- 
ла (21) метода обратной вероятности рассчитьваєт козффициентьт максимального правдоподо- 
бия по сравнению с зталонньми аномалиями (рассчитанньми от теоретических моделей). На- 
дежно вьішявленньми считают аномалиий с вероятностью более 5090. В нашем случаєе вероят- 
ность вьіделения " полезной" аномалии составляет 75-10092, а аномалий-помехи не более 2092. 


Вьгводьг6. На основе анализа магнитньх данньх для парка Ріпеміеуу 30х80 м осуществлен ком- 
плекс преобразований магнитньх аномалий с целью их археологического истолкования. 
Разделение аномалий проведено с помощью аналитического продолжения вверх, а расчет исто- 
чников в зквивалентном слое - с помощью вьтчисления разностного поля от верхней и нижней 
аномальньх границ источников. Дополнительную информацию для качественного анализа дал 
расчет горизонтального градиента магнитной индукции с целью более точной классификации 
локальньїжх аномалий по амплитуде. Зти преобразования можно дополнять расчетами других 
характеристик, информативньх в конкретньїх условиях применения. Результатьт интерпрета- 
ции имеют качественньгй характер, но позволили локализовать места раскопов, рассчитать па- 
раметрью обьектов (намагниченности, глубиньш залегания, контурьш, степень сохранности) с 
помощью решения соответствующих обратньх задач. 
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